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Abstract 

The paper is devoted to the study of a mathematical model for the thermo- 
mechanical evolution of metallic shape memory alloys. The main novelty of our 
approach consists in the fact that we include the possibility for these materials 
to exhibit voids during the phase change process. Indeed, in the engineering 
paper [57j has been recently proved that voids may appear when the mixture is 
produced by the aggregations of powder. Hence, the composition of the mixture 
varies (under either thermal or mechanical actions) in this way: the martensites 
and the austenite transform into one another whereas the voids volume fraction 
evolves. The first goal of this contribution is hence to state a PDE system cap- 
turing all these modelling aspects in order then to establish the well-posedness 
of the associated initial-boundary value problem. 

Key words: shape memory aUoys, mixtures with voids, nonlinear PDEs system, 
initial-boundary value problem, existence of solutions, continuous dependence result 

AMS (MOS) subject classification: 80A22, 74D10, 35A05, 35Q30, 35G25. 

1 Introduction 

Shape memory alloys are mixtures of many martensites variants and of austenite. 
They exhibit an unusual behavior: even if they are permanently deformed, they can 
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totally recover their initial shape just by thermal or mechanical means. There may 
be voids in the mixture, which may appear when the mixture is produced by the 
aggregations of powders, as it has been recently proved in the engineering paper [57] . 
The composition of the mixture varies: the martensites and the austenite transform 
into one another whereas the voids volume fraction evolves. These phase changes 
can be produced either by thermal actions or by mechanical actions. The striking 
properties of shape memory alloys result from interactions between mechanical and 
thermal actions (cf., e.g., [HlllS])- 

We assume that the phases can coexist at each point and we suppose that, 
besides austenite, only two martensitic variants are present. However, this choice 
provides a sufficiently good description of the phenomenon, as we want describe a 
macroscopic predictive theory which can be used for engineering purposes. The phase 
volume fractions, which are state quantities, are subjected to constraints. In particular, 
their sum must be lower than 1, but not necessarily equal to 1, due to the presence of 
voids (cf. also [iO] and [H] where this property is introduced in order to treat solid- 
liquid phase transitions with the possibility of voids and [T] for the corresponding 
numerical results). It is shown that most of the properties of shape memory alloys 
result from careful treatment of those internal constraints (cf., e.g., [Ml^lSn])- AH these 
quoted references are related to a three dimensional model taking the temperature, the 
macroscopic deformation and the volumetric proportion of austenite and martensite 
as state variables. Moreover, let us note that there are not too many mathematical 
models describing phase transitions in which the interactions between different types 
of substances and the possibility of having voids is taken into account: we can quote 
only the two contributions [IHl IH] . 

It is beyond our purposes to give a complete description of the existing literature 
on models for SMA. However, restricting ourselves to the macroscopic description of 
these phenomena, we can refer to the main contributions f36]-[39], [H El [311 [35l [55] and 
[31 [HI [121 UHl UHl [13, [Ml [22] (and references therein) describing full thermomechanical 
models and studying the resulting PDEs from mathematical viewpoint respectively. 
We shall instead focus here on a generalization of the Fremond model for SMA intro- 
duced in [36]-[39]. 

Let us then explain in detail which is the main aim of this contribution, compare 
it with the results already present in the literature, and show the main mathematical 
difficulties encountered. As already mentioned, in this paper we deal with a general- 
ization of the model introduced in [32]- [3^ and later on studied in many contributions 
starting from the pioneering paper [25] , where an existence and uniqueness result has 
been proved for the solution of a simplified problem, where all the nonlinearities in 
the balance of energy are neglected and the momentum balance equation is considered 
in the quasi-stationary form and fourth order terms (related to the second gradient 
theory) are taken into account. In the case when the fourth-order term is omitted an 
existence result dealing with the linearized energy balance equation has been proved 
in [21] , while [22] one can find the proof of the existence of solutions to the linearized 
problem by including an inertial term in the momentum balance. We can report 
also of some results when some or all the nonlinearities are kept in the energy bal- 
ance. The full one- dimensional model is shown to admit a unique solution both in 



Fremond-Rocca 



3 



the quasi-stationary case in [30] and in the case of a hyperbohc momentum equation 
in Existence results have been proved also for the three-dimensional model 

(cf. [23l EHl HH]). Finally, let us mention the uniqueness result for the full quasi-static 
three-dimensional model proved in [T7] and an updated and detailed presentation of 
the Fremond model and related system of equations and conditions, applying to the 
multidimensional case as well, which is provided in [HI [12], [23, Chapter 13], and |39| . 
Let us also point out [TT], [T2] for recent existence and uniqueness results in the three- 
dimensional situation, where the various nonlinear terms arising in the derivation of 
the model are accounted. The large time behavior of solutions is investigated in [26] 
in connection with the convergence to steady-state solutions and in [271 IM] where 
the authors characterize the large time behavior according to the theory of dissipative 
dynamical systems. 

However, all these contributions were dealing with the case in which no voids 
can occur between phases. To model this possibility and to solve rigorously the results 
PDF system is just our aim here. First, in the next Section [21 we derive a model 
taking the possibility of having voids into account, introducing rigorously the pressure 
which has a paramount importance on the mechanical behaviour. In order to do that, 
we follow the ideas of [ID] in which this was done in case of a two-phase transition 
phenomenon. Then, in Section [3l we give a rigorous formulation of initial boundary 
value problem associated with the resulting PDEs and we state our main results: 
existence, uniqueness, and continuous dependence of solutions from the data. The 
proofs are carried over in Section [D and Section [3 The main mathematical difficulties 
are due to the nonlinear and singular coupling between the equations. In particular, in 
order to describe the evolution of the absolute temperature variable, we shall use the 
entropy balance equation (cf. [T3]-[T5] for a complete derivation and motivation of this 
equation). This equation turns out to be singular in {} but the main advantage of using 
it is that once one has proved that a solution component does exists then it turns 
automatically out to be positive and the proof of positivity of the absolute temperature 
is historically one of the main difficulties of these types of problems. The idea here 
is to approximate the nonlinearities with regular functions, to solve the regularized 
system by means of a Banach fixed point argument and then to use compactness and 
lower semicontinuity arguments in order to pass to the limit and obtain a solution of 
the original problem. 

2 The derivation of the model 

In this section we explicitly derive a macroscopic model describing the evolution od 
SMA with the possibility of voids. The model is obtained by properly choosing the 
state quantities, the balance laws and the constitutive relations in agreement with the 
principle of thermodynamics and with experimental evidence. 

2.1 The State Quantities 

We deal only with macroscopic phenomena and macroscopic quantities. To describe 
the deformations of the alloy, the macroscopic small deformation e{u.) , (u being the 
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small displacement) and the temperature are chosen as state quantities. 

The properties of shape memory alloys result from martensite-austenite phase 
changes produced either by thermal actions (as usual) or by mechanical actions. On 
the macroscopic level, some quantities are needed to take those phase changes into 
account. For this purpose, the volume fractions f5i of the martensites and austenite 
are chosen as state quantities. For simplicity, we assume that only two martensites 
exist together with austenite. The volume fractions of the martensites are (3i and (32 ■ 
The volume fraction of austenite is jS^. These volume fractions are not independent: 
they satisfy the following internal constraints 

0<A<1, ^ = 1,2,3 (2.1) 

due to the definition of volume fractions. Since we assume that voids can appear in the 
martensite-austenite mixture, then the /5's must satisfy an other internal constraint 

/3i+(32 + f33<l (2.2) 

the quantity v = I — {Pi + P2 + Ps) being the voids volume fraction. This is the 
case when the alloy is produced by aggregating powders as shown in [57] . In case no 
voids are considered, this sum should be equal to 1 and this considerably simplifies 
the analysis (cf., e.g., PS])- 

We denote by f3 the vector of components Pi (z = 1,2,3) and the set of the 
state quantities is 

E = {eiu),(3,\/(3,'d} 

while the quantities which describe the evolution and the thermal heterogeneity are 

The gradient of V/3 accounts for local interactions of the volume fractions at their 
neighborhood points. 

2.2 The mass balance 

Assuming the same constant density p (the reader can refer to [H] for a model in 
which different densities of the substances are taken into account in a general two- 
phase change phenomenon) and the same velocity U = for each phase, the mass 
balance reads 

piPi +P2 + P3)t + piPl + P2 + P3)divV = . 
Within the small perturbation assumption, this equation gives 

piPi +P2 + P3)t + p{Pl + Pi + /33°)div U = 
where the /3° 's are the initial values of the Pi . For the sake of simplicity, we assume 

Pi + Pi + /?3 = 1 (2.3) 

and have dt{Pi + P2 + Ps) + divU = 0, hence 

dtiPi+P2 + P3) + diYUt = 0. (2.4) 

Mass balance is a relationship between the quantities of 6E , indeed, its effects will be 
included in the cinematic relations (cf. (12.91) in the following subsections). 
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2.3 The equations of motion 

They result from the principle of virtual power involving the power of the internal 
forces, (cf, e.g., [37] ) 

- [ {a: D{V) + B- 5 + H -.Vdjdn 

where V and 5 are virtual velocities, the actual velocities being U and (3^. The 
internal forces are the stress a, the phase change work vector B, and the phase 
change work flux tensor H . The equations of motion are 

pUt = diva + f, = diviy-B + A mfl (2.5) 
an = g, Hn = a on d^l (2.6) 

where p is the density, Ut the acceleration of the alloy which occupies the domain Q , 
with boundary dQ and outward normal vector n . The alloy is loaded by body forces f 
and by surface tractions g , and submitted to body sources of damages A and surfaces 
sources of damage a (in the following we will suppose, for simplicity, A = a = 0). 

2.4 The free energy 

As explained above, a shape memory alloy is considered as a mixture of the martensite 
and austenite phases with volume fractions (3i . The volume free energy of the mixture 
we choose is 

3 

^f = ^{E) = J2 A^i(^) + HE) (2.7) 
1=1 

where the \E'j's are the volume free energies of the i phases and /i is a free energy 
describing interactions between the different phases. We have assumed that inter- 
nal constraints are physical properties, hence, we decide to choose properly the two 
functions describing the material, i.e., the free energy and the pseudopotential of 
dissipation $, in order to take these constraints into account. Since, the pseudopoten- 
tial describes the kinematic properties (i.e., properties which depend on the velocities) 
and the free energy describes the state properties, obviously the internal constraints 
(12. ip and (12. 2p are to be taken into account with the choice of the free energy \E' . 

For this purpose, we assume the ^j's are defined over the whole linear space 
spanned by Pi and the free energy is defined by 

<fiE) = /3i^i(E) + /32^2(^) + (33^3{E) + h{E) . 

We choose the very simple interaction free energy 

h{E) = IcifB) + I \V/3\' 
where Ic is the indicator function of the convex set 



C = {(7i, 72, 73) e < 7, < 1; 71 + 72 + 73 < 1} 



(2.8) 
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Moreover, and by (A;/2)|V/3p we mean the product of two tensors V/3 multiplied by 
the interfacial energy coefEcient {k/2) > 0. The terms Ic{f3) + {k/2)\Vf3\'^ may be 
seen as a mixture or interaction free-energy. 

The only effect of Ic{(3) is to guarantee that the proportions Pi, (32 and (3?^ 
take admissible physical values, i.e. they satisfy constraints (12.11) and (12. 2p (cf. also 
(12.81) ). The interaction free energy term /c(/3) is equal to zero when the mixture is 
physically possible (/3 G C) and to +oo when the mixture is physically impossible 
(/3^C). 

Let us note even if the free energy of the voids phase is , the voids phase has 
physical properties due to the interaction free energy term (z//2)|V/3p which depends 
on the gradient of /3. It is known that this gradient is related to the interfaces 
properties: V/3i , describes properties of the voids-martensites interfaces and 

V/^s describes properties of the voids-austenite interface. In this setting, the voids 
have a role in the phase change and make it different from a phase change without 
voids. The model is simple and schematic but it may be upgraded by introducing 
sophisticated interaction free energy depending on /3 and on V/3 . 

For the volume free energies, we choose 

^i(E) = ^£(u) : Ki : £(u) + ai{^) : e(u) - C^^logd, 

^2{E) = ^eiu) : K2 : ^(u) + a2i^) : £(u) - Ca^^log^?, 

*3(£^) = ^£(u) : Ks : £(u) - -^(^ - ^q) - Cs^^log^, 

where Ki are the volume elastic tensors and Ci the volume heat capacities of the 
phases. Stresses (Ji{'d) depend on temperature ^ and the quantity la is the latent 
heat martensite-austenite volume phase change at temperature i^Q (see Remark 12.11 
below) . 

Remark 2.1. To make the model more realistic, we can introduce two temperatures 
to characterize the transformation: ■(^o , the temperature at the beginning of the trans- 
formation and T?/ the temperature at the end. The interaction free energy is completed 
by hif3) = {la/^o)i^o-^fmf (cf. [6], [55], [56]). 

Because we want to describe the main basic properties of the shape memory 
alloys with voids, we assume that the elastic matrices Ki and the heat capacities Ci 
are the same for all of the phases: 

Ci = C, Ki = K z = 1, 2, 3 . 

Always for the sake of simplicity, we assume that 

(Ti{^) = -(T2{^) = -rW 

where I stands for the identity matrix. Concerning the stress T{i}) , it is known that 
at high temperature the alloy has a classical elastic behaviour. Thus r(?9) = at high 
temperature, and we choose the schematic simple expression 

t{^) = {^- i^,)t, for 7? < r{^) = 0, for 1^ > 



Fremond-Rocca 



7 



with r < and assume the temperature 'dc is greater than . With those assump- 
tions, it results 

ME) = '^'-"^'^^'^ {.(u) : K : e(u)} 
-(A - ■■ e[n) - .ft-^f (i) - »„) - Ct>\ot » + ^ |V/3|' + IdP) . 

2.5 The pseudo-potential of dissipation 

The dissipative forces are defined via a pseudo-potential of dissipation $ introduced 
by J.J. Moreau (it is a convex, positive function with value zero at the origin, [33], 
j49j . |5Uj). As already remarked, the mass balance (12. 4p is a relationship between 
velocities of bE . Thus we take it into account in order to define the pseudo-potential 
and introduce the indicator function Jq of the origin of M as follows 

/o(5i(A+/?2 + /53)+divUt). 

^From experiments, it is known that the behaviour of shape memory alloys depends 
on time, i.e., the behaviour is dissipative. We define a pseudopotential of dissipation 

$(^9, A, V^) = 1 1/3/ + ^1 VAp + A I v^|2 + /o(9*(/?i + /32 + dz) + divu,) (2.9) 

where A > represents the thermal conductivity and c > 0, t; > stand for phase 
change viscosities. 

2.6 The constitutive laws 

The internal forces are split between non-dissipative forces a"'^, B"*^ and iJ"'^ de- 
pending on [E, x, t) and dissipative forces by {a'^, B'^, , — Q'^} depending on bE = 
{£(u)i,/3„V/3i,V7?} and {E,x,t) 

a = a""^ + B = B"'^ + B^ E = H""^ + H'^ 
with the entropy flux vector Q being 

Q = Q''. 

The nondissipative forces are defined with the free energy 

cj'^^E) = -^(E) = (A + P2 + P3)K : e{u) - (A - P2)rW (2.10) 



1 

B"^(E) = T77^(E) = - 

^ ^ d(3^ ' 2 



£(u) : K : e{u) - 2t{^) : e{u) 

K : £(u) + 2r(^9) : e(u) | + B"'^"(E) (2.11) 
K : e{u) - 2J-^^{^ - ^0, 

B"*(E) G dIciP) (2.12) 

i/"^ = kVl3 (2.13) 
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and the dissipative forces are defined with the pseudo-potential of dissipation 

{a^, H^, -Q^} = d<l>{E, 5E) (2.14) 
where the sub differential of $ is with respect to 5E . Relationship (12.141) gives 



a'^ = -pi 


(2.15) 






= 1^ 1 j + 


(2.16) 








(2.17) 




(2.18) 



where p is the pressure in the mixture and it results 

-pe dloiWi +(32 + (3^) + divui) . (2.19) 

The state laws fl2.10l) - fl2.13l) . besides implying that the internal constraints are sat- 
isfied, give also the value of the reactions, during the evolution, to these internal 
constraints. 

Relationships (I2l0l) -( l213l ) and (l2:T5l) -( l2l9D give the constitutive laws 

cr={(3i + (32 + (3^)K : £(u) - ((A - /?2)r(^) + p) I, (2.20) 

{e(u) : K ■ eiu) - 2r(i9) : e(u) - p 1 
e{u) : K : £(u) + 2r{d) : ^(u) - p L ^ndr^j^^ ^ ^^^^ ^2.21) 

B"<^'^(E) G dlcil3) (2.22) 
-p G dIo{dt{(3, +132 + (3^) + divui) (2.23) 
H = kV(3 + vVf3t (2.24) 

-Q(E, 6E) = -Q^E, 5E) = -^W'd . (2.25) 

V 

It can be proved that our choice is such that the internal constraints and the second 
law of thermodynamics are satisfied (cf., e.g., [SHI EZ] and the next Subsection 12. 7p . 

2.7 The entropy balance 

By denoting 

^ = -|^ = C^(l + log^)+/33^ (2.26) 

the entropy balance is 



. -1. , 1 f . /, N , T.d9f3 , rrd 



- + div Q = + - : e{n)t + B'^^ + H"" : V(3, - Q -V^ 

= i?+i|c|Ap + t;|V/3,p + ^|V^p|, inn (2.27) 
-Q ■ n = n, inQ (2.28) 
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because 

due to fl2.19l) . i^Q is the heat flux vector, Ri) is the exterior volume rate of heat that 
is supplied to the alloy, ^^vr is the rate of heat that is supplied by contact action, e{ut) 
is the strain rate. The constitutive laws, within the small perturbation assumption 
and (12.31) . become 

a = K: e{u) - ((A - P^M^) + p) I (2.29) 
-p G diom +132 + f33)t + divU) (2.30) 
-r{^'}):e{u)-p 1 

r(^):5(u)-p Lb«'^^ + c/3, (2.31) 

B""^" G dlcif3) (2.32) 
H = kVl3 + vVf3t (2.33) 

Q(E, 5E) = q,\E, 5E) = -\w'd . (2.34) 

V 

2.8 The set of partial differential equations 

We assume also quasi-static evolution and, using again the small perturbation assump- 
tion, we get the following set of partial differential equations coupling the equations 
of motion (12.51) . the entropy balance (12.271) and constitutive laws (I2.29p -( l2.34p 

div {{K : £(u) - ((/3i - (32)r{'d) + p) I)) + f = (2.35) 

-p G dIo{dt{(3^ +(32 + Ps) + divui) (2.36) 

/-r(^):£(u)-p\ 
c(3t - vA(3t - kA(3 + t{^) : e{u) -p \ + B"'^" = (2.37) 

B"'^'^ G dlc{(3) (2.38) 

+ ^d,(3s - AA log ^ = R. (2.39) 

Ot 

This set is completed by suitable initial conditions and the following boundary condi- 
tions: 

an = g on Si := Ti x [0, T] (2.40) 

u = = on So := To X [0, T] (2.41) 

dnP + dnl3 = on S := dn x [0, T] (2.42) 
dn{\n^) = U onS (2.43) 

where dn is the normal outward derivative to the surface dQ , g is the exterior contact 
force applied to Fi, where (Fq, Fi) is a partition of dQ and Fq, Fi have positive 
measures. 
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2.9 Remarks on the model 

The evolution of a structure made of shape memory alloys, i.e., the computation of 
E{x,t) = {e{u){x,t), (3i{x,t), P2{x,t), i33{x,t),{}{x,t)) depending on the point x of the 
domain Q occupied by the structure and on time t, can be performed by solving nu- 
merically the set of partial differential equations resulting from the equations of motion 
fra . diSD the energy balance flTTTI) . fl^:^ and the constitutive laws (^M>-(^M>, 
completed by convenient initial and boundary conditions (cf., e.g., ^30j, ^J, [61j, 
[53]). The model we have described here is able to account for the different features 
of the shape memory alloys: in particular, their macroscopic, mechanical and thermal 
properties. We have used schematic free energies and schematic pseudopotentials of 
dissipation. 

There are still many possibilities to upgrade the basic choices we have made to 
take into account the practical properties of shape memory alloys. Let us, for instance, 
mention that the pseudopotential of dissipation can be modified in order to describe 
more precisely the hysteretical properties of the materials. There is no difficulty in 
having more than two martensites, for instance, to take care of 24 possible martensites! 
In the same way, it is possible to take into account of the different forms of a single 
martensite variant, as explained in [55] . 

Note that the physical quantities for characterizing an educated shape memory 
alloys are K , C , la, i^^o, ^c, t , the two martensite volume fractions (for the free 
energy) and c, k, \ (for the pseudopotential of dissipation). They are indeed not 
so many in order to have a complete multidimensional model which can be used for 
engineering purposes. 

Other models and results may be found in [H [3, HHl [52], HI] . 

Let us also note the very important role of internal constraints and of the 
reaction B"*' to those internal constraints which are responsible for many properties. 

3 Main results 

In order to give a precise formulation of our problem, let us denote by a bounded, 
convex set in M" (n = 1, 2, 3) with Lipschitz boundary F, by T a positive final time, 
and by Q the space-time cylinder Q x (0,T). Let (ro,ri) be a partition of dQ into 
two measurable sets such that both Fq and Fi have positive surface measure. Finally, 
denote by E := dfl x [0,T], := F^ x [0,T] (j = 0,1) and introduce the Hilbert 
triplet (V, H, V) where 

H:=L\Q) and V := W^'\Q) (3.1) 

and identify, as usual, H (which stands either for the space L'^{Q) or for (L^(fi))^ 
or for (L^(f2))^) with its dual space H' , so that V "—>■ H "—>■ V with dense and 
continuous embeddings. Moreover, we denote by || • \\x the norm in some space X 
and by (■, ■) the duality pairing between V and V and by (■, ■) the scalar product in 
H. 

Set, for simplicity of notation and without any loss of generality 
c=k = la/ A = C = X = v = 1. 
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Then, in order to write the variational formulation of our problem fl2.35f[^7SUl) . 
we need to generalize the relationship B"'''' G die {13) stated in fl2.38l) (cf. also for 
similar generalizations). Hence, we need to introduce the following ingredients 

j : ^ [0, +oo] a proper, convex, lower semicontinuous function 



such that j(0) = and its sub differential (3.2) 

a = d]:W' ^ 2^' (3.3) 

reiy^'°°(M). (3.4) 

Moreover, we consider the associate functionals 

J^(v) = I j(v(x))c/x if V e and j(v) e L\n) (3.5) 

J^(v) = +oo if vGi/ and ^ L\n) (3.6) 

Jy(v) = Jh(v) if V e ^3 (3.7) 

with their subdifferentials (cf. [7, Chap. 11, p. 52]) 

dvyJv : ^ 2(^')' (3.8) 

and (cf. dSl Ex. 2.1.4, p. 21]) 

dnJH : H ^2". (3.9) 



Denote by D{dv,v'Jv) '■= {v G V^^ : dyyJviy) 7^ 0} the domain of dvyJv- Then, 
for X, ^ G H , we have (see, e.g., [T6| Ex. 2.1.3, p. 52]) that 

^ G 9jj Jij(X) if and only if ^ G a(X) a.e. in fl 

and, thanks to (13. 7p and to the definitions of dvyJv and duJu, we have 

dnJuiX) CHn dvyJvW VX G 1/1 (3.10) 

Now we denote by W the following space 

W := {v eV^ : V = on Tq} (3.11) 

endowed with the usual norm. In addition, we introduce on W x W a bilinear sym- 
metric continuous form a(-, ■) defined by 




Note here that (since Fq has positive measure), thanks to Korn's inequality (cf., e.g., 
[20] . [32| p. 110]), there exists a positive constant c such that 

a(v,v) > c||v]|^ VvGlF. (3.12) 
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Next, in order to rewrite the problem (I2.35I - [?!H^ in an abstract framework, let us 
introduce the operators 

B-.V-^V, {Bu,v) = [ Vu-Vv u,veV (3.13) 

A:W-^W', vy'(^u,v)H/ = a(u,v) u,vGiy (3.14) 

H:H^W', w'{'Hu,y)w = I udiyw ueH,weW. (3.15) 

Jn 

Moreover, we make the following assumptions on the data 

uo eW, e L\n), (3.16) 

# > a.e. in fi, wq := log# G W'^'^Q), (3.17) 

/3o = iPlPlPl) e D{dv,v'Jv), (3.18) 

f e W''\0,T;H), g G W''\0,T; {L^iT,)f), (3.19) 

i?G iJ^'^/^(g) ^L^(o,^;L°°(^])), n g l°°(s) □//^'/^'^/^(s). (3.20) 

Then, we introduce the functions 7^ G L2(0,T; 1/') and T G iy^'HO,T, VT') such that 

{n{t),v)= R{t)v+ Il{t)v\Qn veV, for a.e. t G [0,T], (3.21) 
Jq Jan 

w'{0'{t),^)w= I m-^+ I gW-V|an vGl^, for a.e. t G [0,T]. (3.22) 
Now take the function 

7(r) := exp(r) for r G M (3.23) 

and take Jy as in (13.71) . then we are ready to introduce the variational formulation of 
our problem as follows. 

Problem (P). Given > find (u, w, f3i, (^2,(^3) and (^1,^2,^3, p) with the regu- 
larities 

u G L°°(0,T,Vr), divu G iy^'^(0,T;1/) (3.24) 

w G W^'\0, T; H) n L2(0, T; V) H L°°(Q) (3.25) 
/3i,/32,/53eiy^''(0,T;y)nL°°(0,T;l^), f3 e D{dv,v'Jv) a.e. in(0,T) (3.26) 

ei,6,e3eL2(o,r;n, peL2(Q) (3.27) 

satisfying 

^u-7^(p+ (/5i -/52)r(7(«;))) = in VT' a.e. in [0, T] (3.28) 

dtiPi +P2 + Ps) + divut = a.e. in Q (3.29) 

dtw + dtiP-s) + Bw = n inV a.e. in [0, T] (3.30) 

/Bdt(3A /B/3,\ f-ri^iw)) : e(u) - p^ 

dt/3+ Bdtl32 \ + \B(32\+$.= \ r(7(w)) : £(u) -p 

V -(7(^) - ^0) -P 

in (V)^ a.e. in [0,T] (3.31) 

^ = (ei,6,e3)e9y,y.Jy(/3) a.e. in[0,T] (3.32) 
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and such that 

u(0) = uo a.e. in fl (3.33) 

w{0) = Wq a.e. in Q (3.34) 

/3(0) = /3o a.e. inn. (3.35) 

Remark 3.1. Obviously we can take j = Ic in f l3.2p with C as in 02.81) (cf. [T, Ex. 3, 
p. 54]) and recover the problem already stated in fl2.35ti2l39l) as a particular case of our 
more general formulation. Note moreover that we can write down equation (13.311) only 
in {V')^ (and by consequence we need to introduce the notion (13. 8p of sub differential 
in (V' f ) because of the V regularity of {BPi,Bp2,Bp3) in flCTD . The difficult point 
in the proof of our result will be in fact the passage to the limit in the two non-smooth 
nonlinearities in equation (13.311) . By the contrary, we aim to remark that this viscous 
term in (13.311) gives more spatial regularity to dtf3 which furnish more regularity to w 
in (I3.30p (cf. the following Theorem 13. 2p and consequently to -d = 7(w) in (I3.3ip . This 
regularity is needed in order to prove well-posedness for our problem. However, we 
can also notice that, from the mechanical viewpoint, since we have introduced in the 
model the elastic (non-dissipative) local interaction term —A/3 (in (12.360 ). it seems 
also reasonable to include in the model the dissipative local interaction term — A/3j . 

We are now ready to state our main result which is the following global existence 
and uniqueness theorem. 

Theorem 3.2. Let the assumptions on the data A3.^3.22\) hold and let T be a positive 
final time. Then PROBLEM (P) has a unique solution on the whole time interval [0, T] . 

Moreover, in the last Section [5], we will get a proof for the following continuous 
dependence result for Problem (P). 

Theorem 3.3. Let T be a positive final time, (uq, Wq, /3q) (i = 1,2) be two sets 
of initial data satisfying conditions ( l3.J6l - f3T5j) . {7l\ JF*) (i = 1,2) be two data of 
Problem (P) satisfying assumptions I13.21U3.2^) with (P, g\ R\ IT) (i = 1,2) as 
in Let {u\ w\ I3\, (51, Pi) (i = 1,2) be two solutions of Problem (P) 

corresponding to these data. Moreover, besides conditions I13.2[]3.3\] . suppose that the 
following hypothesis 

a e C^'\R^) (3.36) 

holds. Then, there exists a positive constant M , depending on the data of the problem, 
such that the following continuous dependence estimate 

3 

II 1 2||2 I II 1 2||2 I 11/^1 /^2||2 

\\W —W \\L°°(0,t;H)nL^{0,t;V) + 11^ ~^ \\w^-^{0,t;W) \\Pj — Pj \\w^''^iO,t;V) 

< m[h - wlWl + K- ulwi. + J2 ll/3f - PTWl 

+ — -F^||i2(o,t;VF') + 11^^ ~ '^^llL2(o,t;y') j (3.37) 
holds for any t G (0, T) . 
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Remark 3.4. Let us note that in this paper we can treat the difficult couphng be- 
tween the phase-equations fl3.3ip in which appears the temperature {} (= ■y{w)) and 
the entropy balance equation (13.301) in which only the function logi) (= w) plays 
some role, using the L°°-bound on the w-component of solution to PROBLEM (P) 
(cf. (13.251) ). Indeed it was just due to the lack of regularity of solutions that in [13] 
(where there was the same type of coupling without the AdtP-term in (13.311) ) the 
authors did not obtain uniqueness of solutions (cf. also [131 Remark 5.2]). 

However, let us, finally, observe that the main advantage of taking the entropy 
balance equation instead of the internal energy balance equation is that once one has 
solved the problem in some sense and has found the temperature := 7(w) , it is 
automatically positive because it stands in the image of the function 7 (cf. (13.231) ). 
Indeed in many cases it is difficult to deduce this fact only from the internal energy 
balance equation (cf., e.g., [22] in order to see one example of these difficulties). Let 
us note that within the small perturbations assumption the entropy balance and the 
classical heat equation are equivalent in mechanical terms (cf. [T3l fl^ Il5])- 



4 Proof of Theorem 3.2 



The following section is devoted to the proof of Theorem 13.21 First we approximate 
our Problem (P) by a more regular PROBLEM (P^), then (fixed £ > 0) we find 
well-posedness for the approximating problem using a iterated Banach contraction 
fixed-point argument and then we perform some a-priori estimates (independent of 
e) on its solution, which allow us to pass to the limit in Problem (P^) as e \ 0, 
recovering a solution to Problem (P). 



4.1 The approximating problem 

We take a small positive parameter e > and approximate dvyJv in (I3.31f[3l32l) , 
let us take the Lipschitz continuous Yosida-Moreau approximation = df = (j^)' 
(cf. [ini Prop. 2.11, p. 39]) of a and the associated functional JH,eiy) = f^j^iyi^)) dx, 
whose differential {dHJH,e) is the Yosida-Moreau approximation of dnJu (cf. [IS, 
Prop. 2.16, p. 47]). Now we aim to recall some properties of this approximation which 
will be useful in order to pass to the limit as e \ 0. Note that the proof of the 
following lemma is a consequence of [21 Thm. 3.20, p. 289], [21 Thm. 3.62, p. 365], and 
of the Lebesgue theorem of passage to the limit under the sign of integral. 



Lemma 4.1. If dvy'Jv,e -^s the Yosida-Moreau approximation of dyyJv (cf- 
Thm. 2.2, p. 57]), then the following inclusion 



dHJH,e ^ dvyJv,e 



(4.1) 
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holds true. Moreover, the following properties hold true (for e \0) 

Jv,e sup Jy^g in the sense of Mosco (cf Def. 3.17, p. 295]), (4.2) 



e>0 



e>0 £\0 



n Jn 



sup Jy^e = liiii — / hence (4.3) 



Jv,£ Jv In the sense of Mosco and (4.4) 

Vv G V^, Vvg V weakly in Jv(v) < liminf Jv,e{^e)- (4.5) 

Finally, for e \0, it holds 

(u, Jy^e(u)) dyy'Jv 1^ the graph sense ( cf |1[ Def 3.58, p. 360]). (4.6) 

Then, let us call 7"^ the following Lipschitz continuous approximation of the 
function 'j{w) = exp(w7) , i.e. the function 

, I expr if r < l/e . ^ ^, 

7=(r) := < ~ (4.7) 
I (r — l/e) exp(l/£) + exp(l/£:) if r > l/e. 

Moreover let S'^ be the inverse function of 7*^, i.e. 

5^(7^(r))=r Vr G M (4.8) 

and let 7"^ be a primitive of the function 7"^ , i.e. 

7^(r) = 1 + / 7^(s) ds Vr G M. (4.9) 
Jo 

Then the following properties of 7^ hold true (cf. also [TH| Lemma 5.1]). 

Lemma 4.2. There holds 

7'(r) > 7'(r), r(5")'(r) > 1 Vr G M. 

We are ready now to introduce the approximating Problem (P^) as follows. 
Problem (P^). Given i^o > find (u^, ly^, /3f, /3f jp"^) with the regularities 

e L°°{0,T,W), divu' eW^'\0,T;V), p'eL\0,T;H) (4.10) 

w' G W^'\0, T; H) n L2(0, T; 1/) n i:°°(Q) (4.11) 

/3^/3|,/33^Giy^'2(0,T;l^) (4.12) 

satisfying 

^u" - n{f + {(31 - f^2)^{n'{w'))) =T \^W' a.e. in [0, T] (4.13) 

dt{i3\ + + /^g^) + div (u^)t = ep^ a.e. in Q (4.14) 

^tw" + dt^l + Bw' = n in V a.e. in [0, T] (4.15) 

/Bdt(3l\ (E(3l\ /-t{Y{w')) : £(u^) - f 

dt(3'+ BdtPI + B(3I +dHJH,sm = r{Y{w')) : e{u^)-f 
\BdtPlJ \B(3l) \ -{Y{w')-^o)-f 

in (y')^ a.e. in [0, T] (4.16) 
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and such that 

u^(0) = uo a.e. in n (4.17) 

w%0) = Wo a.e. in Q (4.18) 

/3^(0)=/3o a.e. infi. (4.19) 

Concerning this approximating problem, we prove hereafter the following well- 
posedness result. 

Theorem 4.3. Let the assumptions l{3.2[\3.22^) hold true. Let T be a positive final 
time and e > 0. Then the Problem (P*^) has a unique solution in [0, T] . 

Proof. Here we are going first to prove local existence (and uniqueness) in a fi- 
nite time interval [0, f| for some t G [0, T] , then we will extend the solution to 
the whole interval [0,T] proving global existence (and uniqueness) of solution to 
Problem (P^) . Hence, let us take i G [0, T] (we will choose it later) and denote 
by X := (W^'^{0,i;V)Y . Fix for the moment {Pi^ , , ^3^) G A:' in the equations 
(l4.13lHrT^ . then, by well-known results (cf. also fl3.19f[HT^(l ). we find a unique ly*^ = 
l2{Pi\ 02% Ps") e W^'^iO, T- H) n C0([0, T]-V)n ^^(O, T; W^^\n)) solution of ^M- 
By [32], Thm. 6.2, p. 168], it is possible to find a unique = 7i(/3i ,^2,^3 ,w^) G 
L°°(0, T; W) solution of (ICTHmH such that div {u')t G ^^(O, T; H) . 

Moreover, if we take these values of and w'^ in (14.161) . we can find a solution 
iPlPlPI) e X (depending on =: T,i0,' J2" Js, w') and =: T2i0,' J2 
of the equation (14.161) again by standard results. 

In this way, we have defined an operator T : X ^ X such that (/5f , /3|) =: 
^(A ) 5/^3 ) • What we have to do now is to prove that T is a contraction mapping 
on X for a sufficiently small i G [0, T] and moreover, repeating the procedure step 
by step in time (this is possible thanks to the regularities properties of the solution 
listed above), we can prove well-posedness for the Problem (P^) on the whole time 
interval [0, T] and conclude the proof of Theorem 14.31 In order to prove that T is 
contractive, let us proceed by steps and forget of the apices e. 

First step. Let [P^ ,02% P3) e X , = T,{P^\ 02% , = r2{P,\ 02 Js) , and 
(01,02,03) = T{0i,02,03 ) {i = 1,2,3). Then, writing two times dHS]) with 03 
{i = 1,2) instead of dt{03) , making the difference, testing the resulting equation with 
{w^ — w^)t, and integrating on (0,t) with t G [0,T], we get the following inequality 

(4.20) 

for some positive constant Co independent of t. Hence, we get, for all t G [0,T], 

Wiw' - ^^)t||i2(o,;H) + ^Ycomw) < Cill(/53)t - iPl)t\\hio,t;H) (4.21) 

being Ci := CqQ^'^ ■ 

Second step. Let us take -i?* = 7^(w*) and ep* := dt{0l +01 + 01) + div (u*)^ {i = 1, 2) 
and write (14.131) with /3| and l3\ , make the difference, test the resulting equation with 
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(u^ — u^)t, integrate on {0,t) with t G [0,T], and use equation (I4.14p . getting the 
following inequality 

Wit)-^'mW+4p'-pYLH0,t;H) 



3 

< 



- ^ H^') m - Pi) - {f^i - m \^(^p" - p") - - 



,1 „.,2||2 



+ C3texp(2/£)||t^;i - w;1|^o([o,t];y) • (4.22) 

Third step. Write equation fl4.16p for -i?* and u* , make the difference between the two 
equations written for i = 1 and i = 2, and test the resulting vectorial equation by 
the vector {{Pl)t — {Pi)t, {P2)t ~ (/^l)*' WDt ~ {Pi)t)- Summing up the two lines and 
integrating on (0, t) with t G [0, T] , we have (exploiting the Lipschitz continuity of 
dHJH,e (cf. m Pi-oP- 2.6, p. 28])) 



< - 



- p]W - p') 



- / - r(^'))divu^ + (divu^ - divu2)r(^92)) _ (3^)^ 
Jq 

+ [ ((^(^9^) - r(7?2))divui + (divu^ - divu^M^^)) {P^ - 0. 
Jq 

(^1 - - 0i\ + ^ E II (/^D* - (/^l)*lli^(o,*;^.) 
1 ^ 

+ 72 E 11^1 ~ /^j llL2{0,i;H) • 



Moreover, using the definition fl4.2l) of 7"^ and the assumption (13. 4p on r , we get the 
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following inequality 



+ C4texp(2/£)||w^ - It; 



2||2 



ICO([0,t];y) 



Csllu^ - u^||i2(o,T;VF) +texp(2/£:)||ty^ - w 



2||2 



.2 3 



(o,t;y) 



ICO([0,i];/f) 

(4.23) 



Fourth step. Summing up the two inequalities fl4.22p and fl4.23p . two integrals cancels 
out, and using fl4.20p . we get 



2 II 2 

CO([0,t];iy) 



< 



^4^ + ^')Eii(/^l)*-('^l)*iiW) (4-24) 



where does not depend on t. Hence, choosing t sufficiently small (this is our i), 
we recover the contractive property of T . Moreover, applying the Banach fixed point 
theorem to T, we get a unique solution for the Problem (P^) on the time interval 
[0, f| . Due to this estimate it is easy to prove that there exists m E N such that 
is a contraction on X . Hence we have a unique solution on the whole time interval 
[0, T] . This concludes the proof of Theorem 14.31 □ 



4.2 A priori estimates 

In this subsection we perform a-priori estimates on Problem (P^) uniformly in e 
which will lead us pass to the limit as e \ and recover a solution of Problem (P). 
We denote by c all the positive constants (which may also differ from line to line) 
independent of e and depending on the data of the problem. For simplicity, we omit 
the subscript e when it is not necessary. 

First a-priori estimate. Test by Ut, KIM by p, (gUSD by 7^(w) + w, (ITTUD by 
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/3j, sum up the resulting equations and integrate over (0,t) (t G [0,T]). The result is 

Jq \ ^ / JQ JQ Jo 

+^ |V/3(t)r + JnAm) + ^l|u(t)||^ + ^|b||i.(o,,H) = I (l^wo) + Iwl^ 

+^I|v/3oIIh + IwMw + ^ w'{:f, ut)w + ^ (7'(^) + w)) 

- [ dMw [ r{-f%w)) {dwu{dtf3, - dtf32) + - /92)divu,) . 

Jq Jq 

Now, following the line of [131 (5.5)-(5.7), p. 1583], we can deal with the source term 71 
recalling (13.211) and using a well-known compactness inequality (cf. [IHl Theorem 16.4]) 
in this way 

' Rl^w) < f\\R{s)h^^n)h%w{smmn)ds (4.25) 

n7^(^)<ci|n]U^(2)l|(7^(u;))^/2lU2(s) 

Jdn 

<^l|V(7^(t.))^/1?^.(o,,^))3+c||(7^(^))i/li.(o,,^). (4.26) 
Moreover, using assumptions (I3.19P and (I3.22p . we get (integrating by parts in time) 

/ w'{J^,ut)w =- w'{J^t,u)w+ w'{^{t),u{t))w - W/'(J^(0),u(0))i4/ 
Jo Jo 

1 /•* 

<c+-||u(t)||^ + miw'Mw (4.27) 

Now, collecting estimates (14.2511^^271) . using Lemma 112] and (14.141) in order to estimate 
the term containing div and employing assumptions (13.41) on r and (13.1611^1^^ on 
the data, we get the following inequality 

+jHAm) + \\Hmw+e\\p\\i,^,^,.,j,)<c+ [ w'+ [ Mi^+c [ 

^ Jq Jo Jo 

+j [ Ml + cW{^)Y'YLHo,t-,H) + m\w'\Mw 

which - via Gronwall lemma - and choosing e small, leads to the following first estimate 

+4p\\hio,t;H) + \\m\hi^o,fy) + l|V/3(t)||^ + JhMH)) < c. (4.28) 
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Then, applying the standard regularity results for linear parabolic equations to (14.151) . 

we get 

||w^||iyi'2(o,T;H)nL°°(0,T;y)nL2{o,T;W/2.2(f7)) <c. (4.29) 

Now, thanks to assumptions (I3.17P and (13.201) . we may apply [IHl Lemma 4.3] to 
equation (I4.14p . getting the further bound 



||w||loo(q) < c. 
Thanks to (14. 7p we can immediately recover 

||7'HIU-(Q) < c 
and, due to Lemma and (I4.28110T1) . we have 



|V7"HI'<c / rZllHf <c / V5^(7'H)V7'H 

Jq 

< c I iVwP < c 



(4.30) 



(4.31) 



and 



\il'{wM<c 



7^(w) 



<c / |(5^)'(7^H)(7^H),| 



< C / \Wt\ < c. 

Iq 



Hence, from these two inequalities it follows that 

\\l'^iw)\\w^'2{0,T;H)nL2(0,T;V)nL'^{Q) < C. 



(4.32) 



Second a-priori estimate. In order to pass to the limit (as e \ 0) in (I4.16p . we need to 
pass to the limit in p*^ . We will use the following Lions' lemma which is stated in this 
form, e.g., in [HOI Rem. 1.1, p. 17] (its proof is due to |1HI Note (27), p. 320] in case 
of a class domain Q and to [51] when Q is only Lipschitz). For further comments 
on this topic the reader can refer to [lOl Remark 4.1]. 

Lemma 4.4. Let Q be a bounded and Lipschitz set in M.^ and let m be a continuous 
seminorni on H and a norm on the constants. Then there exists a positive constant 
c{fl) (depending only on Q) such that the following inequality 



\\u\\h < c{Q){m{u) + ||Vn||v/} 
holds for all u & H with Vu G V . 



(4.33) 



We want to apply this result in order to find the uniform (in e ) bound on . 
First of all let us note that from comparison in (I4.13p . using also the bound (14.280 on 
u with the assumption (I3.19P on JF, we immediately deduce that 



||V/ 



\l'2{0,T;V') 



< c. 



(4.34) 
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Moreover, always by comparison in fl4.13p . we have that 



p^div V 



< c Vv G ly. 



Following the idea of PU], we can choose G such that 



div fv^) dx 



\^ ■ nds 0. 



(4.35) 



(4.36) 



on 



Note that, since fl is regular (it suffices for to be a Lipschitz domain), we can always 
find a G such that (14.361) is satisfied, because, if we take ^^(x) the ball in cen- 
tered in X G Fi with radius e such that -Be(x) PiFq = and consider the parametriza- 
tion of Fi n B^{x) through the Lipschitz function (xi,X2) ^ {xi,X2,f{xi,X2)) , then 
the normal unit vector associated is 



n 



{d^^ip, -d^^ip, 1) 
JT+WW 



Then, if we take = (0, 0, Q with 

C(y) = 




if \x — y\ < e 
otherwise. 



then G IF and moreover we can show that (I4.36P holds because 

1 1 

> 



where L is the Lipschitz constant of and hence 



div 



c 



■ n (is = . 

ri JrinB,{x) ^/l + iVy^P 



ds^O. 



Take now m in Lemma [4.41 as 



m(v] 



vdiv 



Then, m{v) is a seminorm on H and a norm on the constants because of (14.361) . 
Hence, we can apply Lemma 14.41 to with the choices done above and, thanks to 
flTMlO^ . we get the bound 

\\f\\LHQ)<c. (4.37) 

Finally, by comparison in (I4.16P and using the estimates (I4.28P and (l4.32fH7371) . 
we deduce that also dHJH,s{(3) is bounded in L^(0,T; V) . Then, testing (14.160 with 
and then by using again (I4.28P and (14.3211^7^71) and the monotonicity properties 
of , we get also 

\\f3-\\L^io,T;W'^Hn)) <c ij = 1, 2, 3). (4.38) 

Now it remains only to pass to the limit in (14.1311^71^1) as e \ . This will be the aim 
of the next subsection. 



22 



SMA WITH VOIDS 



4.3 Passage to the limit and uniqueness 

As we have just mentioned, we want to conclude the proof of Theorem 13.21 passing to 
the hmit in the well-posed (cf. Subsection l4.ip Problem (P^) as e \ using the pre- 
vious uniform (in e ) estimates on its solution (cf . Subsection I4.2p and exploiting some 
compactness-monotonicity argument. Let us list before the weak or weak-star conver- 
gence coming directly from the previous estimates and well-known weak-compactness 
results. Note that the following convergences hold only up to a subsequence of e \ 
(let us say \ 0). We denote it again with e only for simplicity of notation. From 
the estimates (l4.28f[OSl) and the property fl4.ll) of dHJH,e, we deduce that 



u weakly star in L°°(0, T; W) (4.39) 

divuf divut weakly in L'^{0,T; H) (4.40) 
weakly in PF1'2(0, T; H) n ^^(0, T; W'^''^{n)) 

and weakly star in (Q) n (0, T; V) (4.41) 

7"(w")^?9 weakly star in iy^'2(0,T;/7)nL°°(0,T;y)nL°°(Q) (4.42) 
^ I3j weakly star in 1^1-2(0, T; V) n L°°(0, T; W^^'^ip)) 

(j = 1,2,3) (4.43) 

^ p weakly in L'^{Q) (4.44) 

weakly in L2(0,T;n (j = l,2,3). (4.45) 

Moreover, employing [58l Cor. 5, p. 86], we get also 

w' strongly in L'^{0, T; V) f] C°([0, T]; H) 

and hence a.e. in Q (4.46) 

r ^ strongly in C0([0, T]; H) (4.47) 

P'^^P, strongly in C°([0,T];F) (j = 1,2,3). (4.48) 



Note that fl4.46fH7^ imply immediately the convergence 

7^(w^) ^ -i? = 7(w) and r(7^(w'^)) — > t{-§) a.e. in Q . 

Moreover, the two convergences f l4.45p and (14. 480 along with the property (14. 6p and 
[21 Thm. 3.66, p. 373] give immediately the identification of the maximal monotone 
graph dvyJv, i-e. 

^ G dv,v'Jvil3) in (V)^ and a.e. in [0,T] 

with ^ = (^1, ,^2, ^3) and (j = 1, 2, 3) are the weak limits defined in (14.450 . All these 
convergences with the identifications made above make us able to pass to the limit 
(as e \ or at least for a subsequence of it) in Problem (P^) finding a solution to 
Problem (P) and concluding in this way the proof of Theorem 13.21 Note that the 
convergences hold for all subsequences Ek of e tending to because of uniqueness of 
solutions. Indeed we may prove it in this way. 
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Consider two solutions of Problem (P) (u*, w\ jSl, jSi^ldi.p') {i = 1,2,3) 
corresponding to the same data. Moreover let us take the mass balance equation 
fl3.29p in the following integrated form 

/3i + /?2 + /?3 + divu = /?i(0)+/?2(0)+/?3(0) + divu(0) a.e. in Q . (4.49) 

Then, integrate equation (13.301) over (0,t) (let us call it 1 * f|3.30p with a little abuse 
of notation) and write down two times equations fl3.28f[Hr^^ . 1 * (13.301) . (I3.3ip with 
(u*, w\ P{, /^35 P*) ) niake the difference between the two equations 1 * (I3.30p , and 
test the result with {w^ — w"^) . Make the difference between the two equations (13.280 . 
test the result with (u^ — u^) . Make the difference between the two equations (I3.3ip . 
written for i = 1 and i = 2, and test the resulting vectorial equation by the vector 
iiPi) ~ iPi)^ iPl) ~ (/^l)' iPl) ~ iPi)) ■ Finally, summing up the three resulting equa- 
tions, integrating over (0,t), with t G [0,T], exploiting the monotonicity of dyyJvy 
using equation (I4.49P in order to get rid of the p-terms, and using the fact that 7, 
defined in (13.230 . is a locally Lipschitz continuous function, -i?* = 7(w*) {i = 1,2), and 
are bounded in L°°{Q) (cf. (13.250 ). we get the following inequality 

- ^')\\lHO,t;W) + - ^YLHO,t;H) + ||1 * («^' - w')it)\\l 

+ E (iK/^l - it) WV + II - P]) \\hiO,t;V)) 




for some positive constant c depending on the data of the problem. Let us notice that 
we have estimated the terms containing the nonlinearity r in (I3.28P and (13.310 on the 
right hand side as follows 

/ {{Pi - PlM^') - {(3', - Pl)r{^')) (div (u^ - u^)) 

+ I ((^(^') - ^(^')) divu^ + r(^92)(divui - divu^)) {{(3l - (3l) - {(3l - (3l)) 
Jq 

1 1 

+c f (1 + iidivui^ + iidivui^) {wi - piwi + wi - mi) ■ 

Jo 

The application of the standard Gronwall lemma together with the regularity (I3.24p 
leads to uniqueness of solutions to Problem (P) and concludes to proof of Theo- 
rem n 



5 Proof of Theorem 13.3 

In this section we give the proof of Theorem 13.31 We will use here the same symbol c 
for some positive constants (depending only on the data of the problem), which may 
also be different from line to line. 
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Then, let us take two sets of data (uq, Wq, (3^), {1Z\ JF*) (i = 1,2) of Prob- 
lem (P) and let (u*, w\ I3\, P^, (31, p^) {i = 1,2) be two solutions of Problem (P) 
corresponding to these data. 

Then, write two times equations (13 .281^3311) with (u*, w\ PI, P\, P\, p*) , make 
the difference between the two equations (13.301) . and test the result with 2{w^ — w"^) . 
Make the difference between the two equations (I3.28p . test the result with 2(u^ — u^)^ . 
Make the difference between the two equations (13.311) . written for i = 1 and i = 2, and 
test the resulting vectorial equation by the vector [{Pl)t — {Pi)t, {P2)t ~ (/^l)*' WDt ~ 

Finally, summing up the three resulting equations, integrating over (0, t) , with 
t G [0, T] , and exploiting the Lipschitz continuity of a (cf. assumption (13.361) ). we get 
the following inequality 

\\n\t) - u'ml + \\w\t) - w\t)\\l +\\w'- wYlHoav) 

+E (iK/^D* - if^])t\\lHo,fy) + i/2iiv(/3j - pfmwi) 

i=i 

< C r (1 + lldivUill^ + ||divU2|P^) 11^1 - 

Jo 

3 

+ E (^11/^1 - f3]\\hmH) + 1/2||(/?1)* - {Pf)t\\Uo,fy)) 
i=i 

+ \\ul -uYw + Wwl- wYh + 1/2 E ||V(/5f - P^Yh 

+c||^i - J'Ylho,T;W') + - -J^'WUo.-yy (5-1) 

Let us notice that we have estimated the terms containing the nonlinearity r in (I3.28P 
and (13.311) on the right hand side using (13.291) and the fact that 7, defined in (13.231) . 
is a locally Lipschitz continuous function, -d^ = ■j{w^) {i = 1,2), and are bounded 
in L°°(g) (cf. (13:251) ). as follows 

/ {{Pi - plM^') - ipl - piM^')) (div (u^ - u^),) 

JQ 

+ / {{ri^') - ri^')) divu^ + r(^2)(divui - divu^)) {{p', - p',)^ - {Pl - p',)^) 

JQ 

{iPl - Pl) (r(^i) - r(^^)) + {iPl - Pl) - iPl - PD) r(^2)) [p] - p])^ 

+ / {{ri^') - ri^')) divu^ + r(^2)(divui - divu^)) {{P^ - p',)^ - {Pl - P^)^) 
JQ 

3 3 

+ c r (1 + lldivul^^ + lldivu^ll^) {\\w' - wYh + - n')r^) . 
Jo 
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Moreover, by adding to both sides in the inequahty fl5.ll) 

3 3 „t 3 



3=1 j=l J=l 



we get the following inequality 

\\vi\t) - VL\t)\\l, + \\w\t) - W\t)\\]j + \\W^ - ^'l|i2(o,i;V) 
3 



j=i 

< c r (1 + lldivuif^ + ||divu2||^) - wYh + - 
Jo 

3 3 



+ — -^^||i2(o^T;W^') + 11^^ ~ '^^llL2(0,t;y') ^ • 

Applying now a standard version of Gronwall's lemma (cf. [TBI, Lemme A. 4, p. 156]), 
we get the desired continuous dependence estimate (13.371) . This concludes the proof 
of Theorem 13.31 □ 
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